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Abstract. We address the question how accurately stellar ages can be determined by stellar evolution 
theory. We select the star with the best observational material available - our Sun. We determine the 
solar age by fitting solar evolution models to a number of observational quantities including several 
obtained from helioseismology, such as photospheric helium abundance or p-mode frequencies. Different 
cases with respect to the number of free parameters and that of the observables to be htted are 
investigated. Age is one of the free parameters determined by the procedure. We find that the neglect 
of hydrogen-helium-diffusion leads to ages deviating by up to 100% from the true, meteoritic solar age. 
Our best models including diffusion yield ages by about 10% too high. The implication for general 
stellar age determination is that a higher accuracy than that can not be expected, even with the 
most up-to-date models. Our results also confirm that diffusion as treated presently in solar models is 
slightly too effective. 
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1. Introduction 

Stellar evolution theory is able to derive the age of stars or stellar systems from a few basic observational data and 
therefore can provide constraints on the evolution of galaxies and the universe itself. Recently, it has become evident 
that improvements in the input physics used for the stellar model calculations lead to considerably reduced ages for 
Pop. II stars as observed in globular clusters (Chaboyer & Kim 1995; Mazzitelli et al. 1995; Salaris et al. 1997), thereby 
bringing into consistency globular cluster ages and cosmic expansion ages. In this context, the question of the accuracy 
of stellar age determinations arises. The errors given for globular cluster ages reflect the observational uncertainties 
only, while those of the stellar models are usually ignored. Chaboyer (1995) and Shi (1995) have investigated the 
influence of various physical assumptions on ages determined from the main sequence turn-off and found possible 
variations of up to order 10%. One of these - the inclusion of non-ideal effects in the equation of state - resulted in 
the reduced cluster ages. 

The influence of the model errors on the derived ages depend, of course, on the dating method used. One therefore 
tries to find ways to minimize the influence of possible sources of error or to determine relative ages, which in general 
are more accurate. For a discussion of several methods for globular clusters, see Stetson et al. (1996). In general, 
the predicted colour and/or visual brightness of the turn-off as a function of age is the decisive quantity, although 
differential age determinators using the colour or brightness difference between the turn-off and some time-independent 
cluster-diagram feature are being used to minimize the influence of systematic offsets. 

A completely independent method to determine stellar ages is used for eclipsing binary systems. One can, for 
example, try to match the positions in the Hertzsprung-Russell- or Colour-Magnitude-Diagram of both components of 
known mass with stellar evolution tracks. A necessary condition for a credible age-determination is then that the two 
error boxes are hit for the same age along both tracks. The classical example for such a successful age determination 
is Al Phe (Andersen et al. 1988) with determined ages that agree better than 10%. Schroder et al. (1997) recently 
have investigated C Aur systems in a different, but related context. In the future, detached eclipsing binary systems 
detected in massive photometric searches (e.g. Kaluzny et al. 1996) will allow to determine globular cluster ages much 
more directly (see Paczyhski 1996a for a review). For that purpose, theoretical luminosity-age-relations for a variety 
of stellar masses and compositions will be necessary (Weiss & Schlattl 1997) that will yield the age (and possibly the 
helium content; cf. Paczyhksi 1996a) of cluster stars of known mass, metallicity and luminosity. Again, the question 
will be: how accurate are such age determinations in terms of the theoretical uncertainties? 

To answer this and other questions concerning the accuracy of stellar ages, Paczyhksi (1996b) suggested to inves¬ 
tigate how accurate we can determine the age of the best-known star - our Sun. In calculating solar models, the solar 
age as determined from meteorites usually is one of the parameters to be fitted by the model. Along with luminosity, 
effective temperature and the metal-to-hydrogen ratio ZjX it determines the model parameters (initial helium content; 
mixing-length parameter); the solar model then provides predictions for the neutrino fluxes, p-mode frequencies and 
other quantities. Since the agreement between models and solar properties determined from helioseismological obser¬ 
vations is very good, this is implicitely taken as evidence that the solar age is appropriate for the models. However, 
there are indications that a slightly different age might lead to even better agreement. Dziembowski et al. (1994) have 
investigated a solar seismic model of only 4 billion years. Although they exclude this model, they also state that a 
much smaller change in age (±5 • 10^ years) cannot be ruled out by helioseismology. Schlattl et al. (1997) find that 
their solar model would agree better with observations, if the solar age would be raised. 

In the present paper, we will turn around the standard solar model approach. A set of observables will be given 
which the models have to match. The solar age will now be a quantity resulting from the best fit to the observables. 
That way, we will be able to determine the solar age from stellar models. The comparison with the meteoritic “true” 
age will give us indications how accurate stellar ages can be determined for the case of the most precisely known stellar 
parameters. For any other star, the accuracy will be lower. It will also define the best set of input physics assumptions 
able to yield the most accurate ages. Since the solar mass and luminosity are among the known parameters, the solar 
case is a stringent test case for age determinations of globular clusters based on detached eclipsing binaries. 

In the next section we will shortly review the numerical code and the input physics assumptions of our solar model 
calculations published in Schlattl et al. (1997). Section 3 will present solar age determinations based on fits to different 
sets of solar quantities, among them interior sound speed and p-mode frequencies. The final section will contain our 
conclusions. 
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2. The solar model calculations 

For the present work, we have used the same solar model code as described in Schlattl et al. (1997). Here, we merely 
recall its basic features and input physics details. 

The opacities are combinations of the OPAL-tables (Rogers & Iglesias 1992; Iglesias & Rogers 1996) with either 
those of Weiss et al. (1990) or Alexander & Fergusson (1994) for the low-temperature regime. For the EOS we either 
use a Saha-type equation or the OPAL equation of state (Rogers et al. 1996). If diffusion is included, the diffusion 
coefficients are calculated according to Thoul et al. (1994). For this work, only hydrogen-helium-diffusion has been 
taken into account. For the layers above r « 20 model atmospheres have been used for the calculations of the p-mode 
frequencies and as the outer boundary condition for the inner model. Since the model atmospheres extend into the 
convective envelope and take into account convection according to the mixing-length description with a parameter 
cTat = 0.5, no constant value for the mixing-length parameter Oint can be used. Rather, a smooth transition has been 
achieved by using the function 

a(T) = /(T) X Oiat + (1 - /(T)) X 
where 

f{T) = (l + e-p(^^)) (1) 

and the parameters Tq and AT being 10 800 resp. 180 K. Oint remains the parameter to be determined by the solar 
model and replaces the global value omlt used in models with standard grey atmospheres. The calculations of Sect. 3.1 
were performed with those standard atmospheres and a constant omlt- 

The composition changes due to nuclear reactions are calculated using a network incorporating the p-p and CNO- 
cycle using the same reaction rates as in Castellani et al. (1994). The calculations are started on the pre-main sequence 
with a homogeneous model powered by gravothermal energy only. Both the spatial and temporal resolution are checked 
for their accuracy and therefore allow a controlled precision of the models. For further details the reader is refered to 
Schlattl et al. (1997). 

Table lists the solar values adopted and the references for the quantities resulting from helioseismological obser¬ 
vations and inversions. 


Table 1. Solar values adopted in this paper and references for quantities inferred from helioseismology. 0 denotes the present 
Sun. Yq is the present surface helium content and Rcz the radius of the convective envelope’s lower boundary, v indicates 
symbolically p-mode frequencies and u the isothermal sound speed. All other symbols have their usual meaning. Errors are la 
errors as given in the original papers. The references are: (1) Grevesse & Noels (1993); (2) Dziembowski et al. (1994); (3) Basu 
et al. (1996); (4) Degl’Innocenti et al. (1997); (5) Christensen-Dalsgaard et al. (1991); (6) Libbrecht et al. (1990); (7) Elsworth 
et al. (1994) 


quantity 

value 

source 

Teff,0 

5777 ±2.5 K 


Le 

(3.844 ± 0.01) • 10®® erg/s 


Rq 

(6.959 ±0.001) ■ 10 “ cm 


iZ/X)e 

0.0245 ± 0.001 

(1) 

Ye 

0.244 

(2) 


0.246 

( 3 ) 


0.238 - 0.259 

( 4 ) 

Rcz 

0.713 ±0.003 Rq 

( 5 ) 

u 


(2) 

V 


(6) 



( 7 ) 

t© 

(4.57 ± 0.02) ■ 10® yrs 
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3. Solar age determinations 

3.1. Solutions of the 3-parameter problem 
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The usual approach to the solar model problem is to ht simultaneously luminosity and radius (or, equivalently, 
effective temperature) with a model having the solar age. In the models, the initial helium content and omlt are the 
two parameters that can be varied and therefore the system has a well-determined solution. In calculations without 
metal diffusion, the metallicity Z is no independent parameter, because it is determined by Yi and the requirement 
that [ZjX^Q has to be matched by the solar model. In calculations taking diffusion into account, Zi is an additional 
parameter and {ZjX )q another solar quantity to be matched. The variations in Z, however, are only of order 1% or 
smaller and not important for our results. We therefore do not count Zi as a further parameter, although in practice 
it has been varied to match {ZjX )q in the calculations presented in this section. 

Accordingly, in this section we solve the following 3-parameter problem: in addition to initial helium content and 
ccmlt (resp. aint) we want to solve for Iq. Therefore, we have to add an additional quantity to be matched. For this, we 
chose either the surface helium content Yq or the depth of the convective envelope i?cz (Tab. p. For both quantities 
we have selected three values repesenting the mean value and the allowed lower and upper limits. 


Table 2. Results of the calculations matching the solar luminosity, effective temperature and a third quantity (second column) - 
either Y© or Rcz/Rq. The third column contains the corresponding other quantity for comparison (Y is the surface He content 
of the final model). Calculations with and without hydrogen/helium-diffusion (forth column) also differ w.r.t. opacities and 
EOS. Standard Eddington grey atmospheres were used, t is the age in 10® yrs as determined. 


Case 

par. 


Diff. 

Qmlt 

F 

t 

A 

Fq 

Rcz / Rq 



A.l 

0.238 

0.693 

n 

2.00 

0.238 

9.00 

A.2 

0.244 

0.700 

n 

1.90 

0.244 

8.09 

A.3 

0.259 

0.720 

n 

1.71 

0.259 

6.16 

A.4 

0.238 

0.722 

y 

1.72 

0.267 

4.34 

A.5 

0.244 

0.727 

y 

1.68 

0.269 

4.08 

A.6 

0.259 

0.743 

y 

1.53 

0.282 

2.64 

B 

Rcz/Rq 

Fs 



B.l 

0.710 

0.253 

n 

1.78 

0.253 

6.86 

B.2 

0.713 

0.254 

n 

1.77 

0.254 

6.71 

B.3 

0.716 

0.257 

n 

1.73 

0.257 

6.38 

B.4 

0.710 

0.233 

y 

1.80 

0.260 

5.15 

B.5 

0.713 

0.235 

y 

1.78 

0.263 

4.81 

B.6 

0.716 

0.240 

y 

1.73 

0.265 

4.48 


Table § summarizes the computational details and the resulting solar ages. For the calculations without 
hydrogen/helium-diffusion we used the Saha-type equation of state (including partial degeneracy) and the Rogers 
& Iglesias (1992) opacities combined with those of Weiss et al. (1990) for the low temperatures. The calculations 
with diffusion made use of the Iglesias & Rogers (1996) plus Alexander & Fergusson (1994) opacities and the OPAL- 
EOS. For all cases listed in Tab. § standard Eddington grey atmospheres were assumed and calculations began on the 
zero-age main-sequence. 

If the helioseismological helium abundance of the photosphere is to be matched (case A), the neglect of diffusion 
leads to much too high ages (up to 97%) due to the low initial helium value, which results in luminosities too low that 
have to be compensated by a higher age. With diffusion, the ages deviate between 5 and 42%. Interestingly, case A.4, 
with the lowest Yq possible is the model providing the best fit, indicating that diffusion might be somewhat too efficient 
in our calculations. This agrees with results by Richard et al. (1997), who found that rotationally induced mixing below 
the convective zone, counteracting diffusion, not only reproduces lithium and beryllium depletion correctly, but also 
improves the agreement in sound speed. 

If the depth of the convective zone is to be matched, again the models without diffusion tend to have ages too high. 
(Note that model B.3 is almost equivalent to A.3.) Models including diffusion determine a solar age that is within 13% 
of the true value. The age of model B.6, which stretches the uncertainties in Rcz/Rq and Yq, deviates by only 2%. 
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It has to be recalled that the models do not include the same input physics as our solar models (Schlattl et al. 
1997) in all respects. However, those including diffusion differ only in the treatment of the atmosphere and the outer 
convective regions. They provide more correct ages, but still deviate by up to 30%. Therefore, we can conclude that in 
the 3-parameter case, given the observational uncertainty in the third parameter, the solar age cannot be determined 
as accurate as 10% for most choices. The exceptions are cases B.4, B.5 and B.6, where Rcz/Rq is matched but Y© 
ignored. 


3.2. Fitting 4 quantities with 3 parameters 

In the previous section it became obvious that the depth of the solar convective zone and the surface helium content 
cannot be fitted simultaneously by one model even when using present state-of-the-art input physics. Consequently, 
we are now going to determine the minimum x^-fit to all four solar quanitites. is defined as 


4 4 

= = ^ 


{A, - B,f 


i=l 


i=l 


( 2 ) 


where Ai resp. Bi denote the observed and model value of effective temperature, luminosity, surface helium content 
and depth of convective zone, and ai the observational errors (see Tab. 




CX) 

CT5 



Fig. 1. Fitting four observational quantities: x = 32 surface and projections onto the t — aint and Yi — aint planes in the 
parameter space 


The models for this set of calculations include, with the exception of metal diffusion, all the physics of the best 
models of Schlattl et al. (1997), i.e. the special treatment of the convective gradient (Eq.(l)) and the inclusion of 
model atmospheres. They also start with the pre-main sequence evolution. Otherwise they agree with those of set B. 

Contrary to the calculations of the previous section, for the present problem we have to scan the 3-dimensional 
parameter space. We did this with the following steps: AYi = 0.001, Aoint = 0.1, At = 2 ■ 10® yrs. For Oint = 5.5, we 
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also used a finer t-grid of 5 • 10^ yrs (see below). As discussed in the previous section we now kept Zi constant at the 
value needed for {ZjX^i = {ZIX)q in order to save computing time. In the investigation by Richard et al. (1997) the 
influence of varying Z on Yi and the resulting solar model can be seen to be very small. We therefore do not expect 
that our results would change if we had not made this simplification. 


Table 3. Models for selected values of Oint with the lowest of the problem fitting four observational quantities. Oint is the 
interior value of the variable mixing-length parameter (see eq. 1). For more details of the models, see the text. Ages are in 10® 
yrs. 


Case 

Q^int 

Yi 

t 

L/Lq 

Teff [K] 

Rcz/Rq 



C.l 

4.5 

0.263 

4.45 

1.0051 

5780.0 

0.7211 

0.2359 

14.84 

C.2 

4.7 

0.261 

4.65 

1.0030 

5777.0 

0.7186 

0.2336 

7.36 

C.3 

5.0 

0.259 

4.85 

1.0010 

5778.0 

0.7157 

0.2316 

3.91 

C.4 

5.1 

0.259 

4.85 

1.0014 

5782.4 

0.7152 

0.2317 

8.27 

C.5 

5.3 

0.257 

5.05 

0.9988 

5777.1 

0.7138 

0.2294 

3.46 

C.6 

5.5 

0.255 

5.25 

0.9959 

5771.3 

0.7113 

0.2271 

11.58 

C.6’ 

5.5 

0.256 

5.15 

0.9977 

5771.2 

0.7119 

0.2284 

4.51 

C.7 

5.7 

0.255 

5.25 

0.9965 

5778.7 

0.7108 

0.2273 

6.41 

C.8 

5.9 

0.254 

5.45 

1.0035 

5781.9 

0.7097 

0.2258 

12.18 

C.9 

6.1 

0.252 

5.65 

1.0001 

5774.5 

0.7070 

0.2235 

10.73 


Figure shows y^-surfaces and projected contour lines for our models. Apparently, there is no clearly defined sharp 
minimum which would allow a unique determination of the solar age from fitting the four quantities. Rather, there 
exists a long-stretched region of equally well fitting sets of parameters (including age). The separated minima visible 
in the projected contour plots are an artifact resulting from the grid in parameter space and the plot interpolation. 
In reality, the = 8 contour lines are connected. We have verified this by using a higher resolution in age for the 
case of ttint = 5.5 (cf. Figs. || and ^). Our grid of models in the aint and Yi parameter space is rather coarse due to 
the otherwise inhibitive number of solar evolution calculations. In Tab. ^ we list those models with the lowest y^ for 
several values of Oint- None of them is fitting all quantities very well (in particular, Yq), as is indicated by the y^-values 
which should be close to 1 for a good model. This is a confirmation of the fact that at the presently given accuracy 
of observations the theoretical models still contain significant deficits. Model C.5 is the best model of the real sun. Its 
age is within 10% of the presumed solar age. C.3 is within 6% of the true solar age at a slightly higher y^. Given the 
uncertainty of « 4% in the age determinations due to the time resolution we use, the true minimum could well be 
within 2-3 % of the solar age. The effect of a finer resolution in At is demonstrated by model C.6’. Relative to C.6 the 
minimum is now located more accurately (y^ = 4.52). The age of the true minimum is between 5.1 and 5.2 Gyr. It is 
evident that all models in Tab. ^ have a surface helium content too low and consequently a rather low initial one, too. 


3.3. Age determination by best-fit models to all observational data 

Up to now we have used only two quantities derived from helioseismology. We are now going to include results about 
p-mode frequencies and sound speed. A simple inclusion of many frequencies or of the sound speed at many radial 
points as additional degrees of freedom would dominate the standard y^-formula completely. In addition, since neither 
the frequencies nor the sound speed are statistically independent from i?cz or Yq, the strict interpretation of y^ is no 
longer valid. We therefore modify Eq.(^ to include two additional terms and interprete the resulting quantity as a 
function of merit, defined as 

= (3) 

where the first four terms are defined as in Eq.(^. The first additional term 

fn(uo(x) — Uth(x))^dx 

W5 = Wn= - 

fo crl(x)dx 


( 4 ) 
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describes the deviation of the model sound speed Uth as a function of fractional radius x = t/Rq from the helioseis- 
mologically inferred one (uo) weighted by the observational error au- The other term 


W6 = W, 




( 5 ) 


similarly is used to calculate the model goodness of the 15 frequencies iy(0,n) for a representative mode. We have 
selected the I = 0 mode because it traces the whole solar interior. The choice of the 15 radial modes was dictated by 
the availability in the original papers. 

Since the observational errors on u and v are extremely small, even the mean deviations Wu and W^, still would 
dominate the function of merit. To allow the global quantities to retain their influence, we normalize the individual 
terms in Eq.(^) by their respective minimum values. 



Fig. 2. Quality of model fits using all solar quantities, illustrated by contour lines of the logarithm of the function-of-merit x^, 
where all contributing terms have been normalized by their individual minimum value. The four panels are (from upper left to 
lower right) for Oint = 5.1, 5.3, 5.5, 5.7. Only parts of the parameter space has been investigated 


For evaluating the same solar evolution calculations as in the last section were used. The result of this age 
determination does not differ significantly from that of case C. In Fig. ^ we show contour lines of logx^, the function- 
of-merit as defined by Eqs.(^|-^, where the individual terms have been normalized by their respective minimum value. 
Evidently, there are only very small regions, where < 10, for all four cases of aint shown. For aint smaller than 
4.7 and larger than 6.1 they vanish completely. The 3-dimensional isocontour resembles that of Fig. ^ very much. 
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5.0 5.2 5.4 5.6 5.8 


Fig. 3. As Fig. 2, but for the case of hxed Yi = 0.257. The minimum at the center corresponds to model C.5 



Fig. 4. Contour-lines of log Wu illustrating the quality of the models with respect to sound speed (fixed Oint = 5.5; fine 
t-resolution). Note that Wu has not been normalized to its global minimum here 


therefore. The effect of a finer resolution in age {At = 5 ■ 10^ yrs) is seen in the lower left panel (aint = 5.5). Here, the 
contours for logx^ < 3 are still elongated, while in all other panels they form isolated wells. We have compared the 
shown aint = 5.5 contour lines with those for the same t-resolution as used for the other three panels {At = 2 • 10® 
yrs) and found the same round shape of the logy^ < 3 contours. The minimum for this case deviates slightly from 
that of model C.6’, which was calculated with the same resolution. It is now at 5.2 Gyr but the same Yi = 0.256. The 
global minimum is reached (see Fig. ^ upper right panel) for aint = 5.3 at Yi = 0.257 and age t© = 5.05 ■ 10® as for 
model C.5. For comparison, a cut at = 0.257 is shown in Fig. The y® minimum again corresponds to model C.5. 

The contribution of Wu and to y® is illustrated in Figs. ^ and Fig. ^ demonstrates that the sound speed 
y®-contours Wu outline a long-stretched valley enclosing the parameters of the best models without defining a clear 
minimum. Rather, the function-of-merit appears to be almost degenerate. Furthermore, the region of smallest y® 
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is shifted towards higher ages and helium content with respect to the location of the best model. The shown case 
(cTint = 5.5) is a typical one. u cannot be used for age determinations. 

Similarly, Fig. ^ shows how the models fit the I = 0 p-mode frequencies. The minima, though not very low, are 
very sharply defined, and the models providing the best fits for the complete problem also have a comparably low 
W,^. Model C.6, e.g., lies within the very narrow \ogW^ = 2 contour of the upper panel, which was plotted using the 
high t-resolution. The minimum within this countour is close to the position of model C.6, but at a slightly larger 
age. The position of model C.6’, which does not include the ht to the p-mode frequencies, but was calculated with 
the same high t-resolution, is outside the log Wi, = 2 contour. Therefore, at a resolution of 5 ■ 10^ years, the minima 
of models C do no longer coincide with those of the p-mode fits. The complete fit (see Fig. ||) incidentally is at the 
position of model C.6. Since in the complete ht we applied an - arbitrary - weighing of the contributing terms to 
by the normalization, the difference between these various minima gives a measure for a systematic age determination 
uncertainty of the order of 2%. 

In the lower panel, Model C.6 can clearly be identihed as the central minimum. The p-modes therefore add very 
important information in selecting the best-ht models and determining the solar age. The also dominates 

and the contours of Fig. |^. We cannot exclude that for an even higher resolution in aint or Yi the mimima deviate 
from those of the x^ or of x^ of models C, similar as in the demonstrated case of a higher t-resolution. 






Fig. 5. Contour-lines of the unnormalized log Wv illustrating the quality of the models with respect to the frequencies of the 
I = 0-mode. Upper panel: for fixed parameter Omt = 5.5 (fine t-resolution); lower panel: fixed Yi = 0.255 













A. Weiss & H. Schlattl: The age of the most nearby star 


9 


4. Discussion 

If the solar age were to be determined from stellar models fitting observed solar quantities, we have shown that a 
priori an uncertainty of a factor of two would exist. The reason for this is the fact that apart from the global solar 
values - radius, effective temperature and {ZjX )q - the third quantity necessary to restrict mixing length parameter, 
initial helium content and age, has to be one derived from helioseismology. In our calculations, this was either the 
depth of the solar convective zone or the present photospheric helium abundance. Both quantities depend crucially on 
the effect of particle diffusion, though does less so (Tab. ||). If this is not recognized, and diffusion not taken into 
account into the calculations, derived ages are by far too high (model sequences A & B without diffusion). If diffusion 
is taken into account, the predicted ages can be as accurate as 10 %, if i?cz is the third parameter, but are too low by 
up to 40% for the allowed range of Yq. The latter fact indicates that either diffusion is too effective in the models or 
that 10 is rather at the lower boundary of the present determinations. 

Including both of these helioseismological quantities in fitting the present sun, x^-fits using “best-physics-models” 
yield ages that agree to 10% or better with the solar age. However, none of the models is a good model for the Sun 
with x^ ~ 1 there is no clear global x^-minimum fit. If we generously accept all models of Tab. I with x^ ^ 10 
(approx. 3 (t), the solar age could range between 4.65 and 5.65 Gyr, being slightly too high. Therefore, a systematic 
error in the models is still present. 

Adding sound speed and p-mode frequencies to the fits, the best-fitting models remain nearly the same as in 
the previous case, if we normalize the individual contributions to the function of merit. This is necessary, since all 
solar models do not predict p-mode frequencies and sound speed throughout the Sun with sufficient accuracy for the 
extremely low errors in the observations. Since the various quantities used in this fit are no longer independent of each 
other, the function of merit, though still denoted by x^j no longer has any statistical meaning. While sound speed does 
not provide helpful information for age determinations, the p-mode frequencies are extremely sensitive to deviations 
of the models from the observed Sun. In fact, though all models are poor fits for the p-modes, the function-of-merit 
minima are extremely sharp and coincide with the parameters also yielding the best fits for the previous case of 4 
solar parameters, as long as we use age-steps of 2 ■ 10® yrs. At a four times higher resolution, the various x^-mininia 
do no longer coincide, but the location of the “true” minimum remains subject to the exact definition of x^- p-mode 
frequencies alone might suffice for finding best-fit parameters, including the solar age. 


Table 4. x^-niinima of the problem fitting four observational quantities for Oint = 5.3 and two values of Y — i. In these 
calculations metal diffusion has been taken into account and the high resolution in time has been applied. For comparison, 
model C.5 is repeated 


Case 

Q;int 

Yi 

t 

L/Lq 

Teff [K] 

Rcz/ Rq 

T, 


D.l 

5.3 

0.257 

5.00 

0.9970 

5779.8 

0.7182 

0.2295 

8.40 

D.2 

5.3 

0.256 

5.15 

0.9994 

5774.5 

0.7180 

0.2277 

9.42 

C.5 

5.3 

0.257 

5.05 

0.9988 

5777.1 

0.7138 

0.2294 

3.46 


We have for one special case investigated the effect of metal diffusion (aint = 5.3; At = 5 • 10^; cf. model C.5). 
Table ^ lists the minima for two different choices of Yi. Compared to C.5, model D.l (same Yi) has a minimum of 
higher x^ at a slightly lower age. A better fit is found for a lower Yi but at a higher age (D.2). At least for this value of 
Oint) we conclude that the consideration of metal diffusion would lead to a x^-minimum at a much higher age. Whether 
the global x^-minimum would also be shifted towards an higher age, can only be decided after the (time-consuming) 
calculations for all parameters have been done. From the fact that the depth of the convective zone has increased for 
models D, we infer that only the C-models with a low value for Rcz/R q could give the global x^-minimum, if metal 
diffusion were included. Those, however, already have ages too high. It might be that the inclusion of “better” physics 
in this case leads to a larger error in the age determination. 

The implications for stellar age determinations are not straightforward. Evidently, the physical input used today 
for the standard solar model (including diffusion) provides the most accurate age determination. It therefore should be 
used in stellar models in general. If direct measurements of stellar parameters (mass, luminosity, effective temperature) 
are available, the derived ages can be expected to be comparably accurate as in the solar case provided the observations 
are as accurate as for the Sun. 

The inclusion of diffusion is the most important aspect of solar age determinations performed as in the present 
paper. If it is neglected, the derived ages are unacceptable. This is due to our knowledge about the helium content 
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or depth of the solar convective zone. For other stars, this information will not be available directly. Evolutionary 
ages will then depend solely on the assumed initial helium content, no matter whether diffusion is taken into account 
or not. Paczyhski (1996a) argues that detached eclipsing binaries will provide the possibility to determine age and 
(initial) helium content at the same time. In this case, diffusion most likely has to be included, since the central helium 
content is changed by diffusion and this influences the main-sequence lifetime. Otherwise, an additional uncertainty 
of about 5% (Chaboyer 1995) is added. Note that for old systems, the relative error due to a wrong helium content is 
smaller than in the solar case, because it leads to a luminosity change roughly constant for the whole main sequence 
evolution. Absolute errors could be of up to a few Gyr. 

For age determination methods based on differential quantities like the brightness difference between turn-off and 
horizontal branch the situation is more complicated. If the defects of the stellar evolution calculations, which lead to 
the solar age errors, affect the different stellar evolution stages systematically, the ages based on differential methods 
could be more accurate than the results found here indicate. (There are, however, many other, probably more severe 
sources of errors, like the conversion from theoretical to observational brightness.) This, however, has to be investigated 
in detail (Castellani et al. 1997). 

To summarize, we have demonstrated that even with the best observational data available, present stellar evolution 
calculations cannot be expected to yield ages with errors less than 10%. Hydrogen/helium diffusion is an absolutely 
necessary ingredient to reach this level of accuracy. We have not investigated, whether the inclusion of metal-diffusion 
would lead to a more accurate solar age, because in most stellar evolution calculations this most likely will not be 
included due to computational limitations. 
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